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Abstract 

We consider the impact of cosmological B — L constraints on supersymmetric standard models with 
bilinear breaking of R-parity. In order to avoid erasing any primordial baryon or lepton asymmetry 
above the electroweak scale, B — L violation for at least one generation should be sufficiently small. 
Working in the context of models with non-universal soft supersymmetry breaking masses, we show 
how the above cosmological constraint can be satisfied while simultaneously providing a neutrino 
mass matrix required by current data. 
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In recent years the increasingly strong evidence for neutrino oscillations from various 



experiments 



l| has led to the active study of R-parity violating extensions of the minimal 



n 



supersymmetric standard model (MSSM) [2]. Such models maintain the particle spectrum 
of the MSSM but contain renormalizable lepton flavour violating couplings. The observed 
neutrino oscillations and mass differences can be accommodated with such couplings 0,0], 
and so these models provide a conceivable alternative to seesaw 

mass generation. In contrast to the R-parity conserving MSSM, the lightest supersymmetric 
particle is unstable and decays in the detector with branching ratios which are correlated with 
the neutrino mixing This provides a robust, experimentally accessible test of the model 
at the Large Hadron Collider and/or a e~^e~ Linear Collider Q]. Analogous confirmatory 
signatures are less readily found for the elegant seesaw mechanism Bilinear R-parity 
violation (BRpV) is the minimal extension of the MSSM with R-parity violating terms 
0, ^, 11 1 . The minimal supergravity version of BRpV (i.e. imposing universal 
soft supersymmetry breaking masses at an ultraviolet scale) can easily accommodate the 
atmospheric neutrino oscillation data. However, in order to provide the currently favoured 
large mixin g an gle solution for the solar neutrino anomaly, this universality condition must 



be relaxed 



13l Il4l Il5l |. Another option for obtaining a realistic neutrino mass matrix is 



to allow both bilinear and trilinear couplings while keeping the universality condition of 
the soft supersymmetry breaking masses. The minimal model of trilinear R-parity violation 
(TRpV) assumes the dominance of the third generation trilinear couplings and thus contains 
five free parameters of lepton number violation to fit all the neutrino data successfully 0]. 
The theoretical background on massive Majorana neutrinos and lepton violating mixing 



matrices describing neutrino oscillations can be found in 



17[. The atmospheric neutrino 



data is exp lained by oscillations z/^ < — > Ur, and a global analysis gives the following 3a 



ranges 



1^ 



0.3 < sin^ ^atm < 0.7 
1.2 X 10"^ eV^ < Aml^^^ < 4.8 x 10"^ eV^ 



(1) 



with maximal mixing sin^ 9^ 



0.5 and Am^^ 



2.5 X 10"^ eV^ as the best fit point. 



Similarly, the solar neutrino data is explained by Uf, oscillation into a mixture of z/^ and z/^. 
Global analysis suggests a large mixing angle, although not maximal, and a much smaller 
mass squared difference. The allowed region for Am^^j previous to KAMLAND results ^| is 



now split into two sub-regions. At 3cr we have [18|, [1 

0.29 < tan^ ^^oi < 0.86 

5.1 X 10"^ eV^ < Am^^i < 9.7 x 10"^ eV^ 

1.2 X 10-^ eV^ < Am^^i < 1.9 x 10"^ eV^ (2) 

with tan^ 6'soi = 0.46 and Am^Qi = 6.9 x 10"^ eV^ as the best fit point. 

In connection with neutrino physics, there appears an important cosmological consid- 
eration. As is well known, the seesaw mechanism provides a natural way to generate the 
baryon asymmetry of the universe through the out-of-equilibrium decay of a heavy right- 
handed neutrino j^]. Being a new physics model just around TeV scale, the R-parity 
violating MSSM can hardly accommodate such a mechanism of baryogenesis. However, in 
the supersymmetric model, the so-called Affleck-Dine mechanism can successfully work to 
generate the required amount of the baryon asymmetry in the flat direction along, e.g., LH^ 



21[. It is notable that such a property is unaltered even with the presence of R-parity 
violating terms which must be very small to generate tiny neutrino masses. 

It is known that lepton number violating couplings have important consequences for 
baryogenesis since together with B + L violating sphaleron processes they are capable of 



erasing any pre-existing baryon/lepton asymmetry in the universe |22l 2^|2^- The purpose 
of this paper is to explicitly check if such cosmological constraints on the lepton violating 
couplings can be satisfled in BRpV while simultaneously accommodating the form of the 
neutrino mass matrix indicated by the atmospheric, solar and reactor neutrino experiments. 
A previous analysis j25| derived the cosmological bounds for BRpV but their effect on the 
neutrino mass matrix was not covered. Given the wealth of new data which has become 
available since appeared, we develop their analysis and apply the bounds to the currently 
favoured bimaximal mixing form of the neutrino mass matrix. 

we note .Ka. ou. .nve..ation . not relevant ,t e.ee.owea. ba„. Q . operative, 
in which case the produced baryon asymmetry cannot be erased solely by R-parity violating 
processes. For our purposes we assume that a. B — L asymmetry was generated primordially 
by some means at a high energy scale, and our intention is its preservation at all energies 
down to the electroweak scale when the sphalerons flnally fall out of equilibrium. 

We briefly summarize the mechanism of neutrino mass and mixing generation by R-parity 
violating couplings, both bilinear and trilinear. The R-parity violating MSSM predicts a hi- 



erarchical neutrino mass spectrum. The atmospheric mass scale corresponds approximately 
to the heaviest neutrino mass, m^, and it is generated at tree level via a low energy see-saw 
mechanism due to the mixing of the neutrinos with the neutralinos. On the other hand, 
the solar mass scale, corresponds approximately to the second heaviest neutrino, m2, and 
is generated at the one loop level . The atmospheric neutrino mixing is also predicted by 
tree level physics, and depends in a simple way on sneutrino vacuum expectation values 
expressed in the basis where the bilinear parameters are removed from the superpotential. 
On the other hand, the solar neutrino mixing angle is again predicted by one-loop physics 
which is mainly determined either by the trilinear couplings in the superpotential or by the 
bilinear parameters in the scalar potential. 

Let us remark, however, that we cannot exclude the possibility of the loop mass domi- 
nating over the tree mass, which may have an interesting implication to baryogenesis as will 
be discussed later. 



The well-known baryogenesis constraint 



22 



23[ can be easily applied to the TRpV model 



with the universality to exclude this possibility. To see this, let us consider the following 
trilinear R-parity violating couplings in the superpotential; 

W = \,uUL,El + K.kUQPl (3) 

which generates a neutrino mass at one-loop level as follows; 

f^fioop _ 3 '^^33'^i33 rnl{Ab + /itan/j) ^ ^m^'i?. ml{Ar + /i tan/3) 

Stt^ m? — m3 rrir Svr^ m% — m? m? 

bi 02 02 ^1 ^2 T2 

Note that we have picked up A and Aj33 which give the largest contribution to the neutrino 
masses when all the trilinear couplings are of similar magnitude. Then, requiring the above 



one-loop mass Q gives rise to the solar neutrino mass scale, m2 ~ y Am'^^^ ^ 8 x 10 ^ eV, 
we obtain 

a;.,awV3.5x:o-(3^)"^(£^)"'^ (5) 

taking m = Ab + fitan/S = "^^i 2 = + yUtan/3 = mf^^. Now, the problem is that such 
a large coupling makes lepton number violating interactions very active when the B + L 
violating sphaleron interaction is also in thermal equilibrium, so together they erase the 
baryon asymmetry before the electroweak phase transition. Indeed, the interaction in Eq. Q 



3 



gives the decay width for lepton number violating one-to-two body decays, 

192C(3) T ^ ^ 

assuming T ^ m. The out-of-equihbrium condition, < H = LGQ^/g^JfT'^/mpi, gives 



for (7e// = 915/4. This is for T » m. An improved resuh which does not make this 
assumption was presented in 2^ and shows that the T/m dependence of Eq. © is very 
mild. One sees a big contradiction between (0) and ((Tj). As indicated in Eq. (0), one needs 
the trilinear couplings of A233333 ~ Ai33^233 ~ C(10~^) to accommodate the required bi- 
large mixing of the atmospheric and solar neutrino oscillation j^. Thus, the baryogenesis 
constraint rules out a purely TRpV explanation of the observed neutrino data. 

The situation may be different if the neutrino masses are generated purely by bilinear R- 
parity violating couplings with non-universal soft masses, in which case the non-universality 
can give much freedom. Forbidding the lepton number violating trilinear couplings in the 
superpotential in Eq. (jHI), the BRpV model allows the following dimension- two terms in the 
superpotential and in the soft supersymmetry breaking scalar potential: 

Vsoft = lJ.{eiBiLiH2 + BH^H^) + ml^jj.LiHl + h.c. (8) 

Here we have used the same notation for the superfields and their scalar components. A key 
point to notice is that without the electroweak symmetry breaking, the S't/(4) rotation in 
the 'superfields', Lj and Hi] 

Li ^ Li + tiHi and Hi ^ Hi - etU (9) 

which gets rid of the term (valid up to 0(ej)) leaves invariant the gauge interactions and 
thus its effect is only to generate the effective couplings as in Eq. Q with 

A -33 = eihb and A^ss = . (10) 

Under the SU{A) rotation Q, the scalar potential in (jH)) becomes 

Vsoft = lJ^iBHiH2 - tiABiUH^) + {m%H, " eiAml)UH\ + h.c. (11) 



where ABj = B — Bi and Am^ = — m\,. Eq. (jlip shows that the additional lepton 
number violating mixing mass terms for the 'scalar fields' Li and Hi^2 (in the basis of vanish- 
ing ej) arise in the presence of the non-universal soft supersymmetry breaking parameters. 
Diagonalizing away such mixing mass terms can be made by the following rotation among 
the scalar fields Lj, Hi and H'2 = iT2H2'- 

Li ^ Li — EiiHi — ej2-f^2 
Hi ^ Hi + EiiLi 

H'^^H'^ + e.^U (12) 

where the variables Sn and £^2 are determined as 

{m\j^ + /x^ - m|J(e^Am^ - mj^j^J - ej^'^BABj 



[m 



^ ^ (rriH, + /i - mjJe./xAE, - fiB{t.Ami - mj^^^J 
("^Hi + /^^ - "^i, ) ("^lf2 + /^^ - "^i J - 
As will be discussed later, it is useful to rewrite £41,22 in terms of the variables and rji 
defined by 

, _ {^^) , _ , ^ AS, 



(i/l) ' " - - ' - 5 

where (i/j) and (Hi) are the vacuum expectation values of the sneutrino and Higgs boson 
generated after the electroweak symmetry breaking. Using the minimization condition of 
the Higgs and sneutrino fields, we obtain 

_ _^ m\sl{m% - M|c2/j) 

tp m\{m\ - m\) - (m| - m\sf)Mlc2(i 

where t/3 = tan/5 = {H2)/{Hi). The variables £ji,i2 control the size of lepton number 
violating interactions which now arise due to the misalignment between the scalars, Li and 
Hi^2, and fermions, Li and ifi,2- That is, the rotation ()12|) gives rise to the following lepton 
number violating vertices: 

^eff ~ hrSixLiL^E^ + hhEiiLiQ^D^ + htei2L[Q^Ul 
+ ^-^[HlUB + L\HiB] + ^-^Uh'2B 

+ ^-^[HyLiX' + L^HiX''] + ^LyiL^X- + h.c. (15) 
V 2 V 2 



where L'- = iT2Ll, are Pauli matrices and A'' represent the Sll{2) gauginos. Applying the 
constraint ^ to the couphngs in Eqs. (|TO|l and (|T5|l . we get 



ei < 1.2 X IQ-^c^ 



m 



300 GeV 



ea < 3xl0-'( V^' (16) 

V300 GeV/ ^ ^ 



V300 GeV. 

where m is the smallest mass of the sfermions involved in the A^33 term; Li, Qs and D^, 
is a gaugino mass involved in the process x ~^ LiHi and the last equation comes from the 
process Lj Q^U^. 

The sizes of certain bilinear parameters are determined to generate realistic neutrino 
masses and mixing in our bilinear model. First of all, upon electroweak symmetry breaking, 
the Higgs and sneutrino acquire vacuum expectation values and generate a tree-level neutrino 
mass matrix 

= (17) 

where Fat = MiM2/(c^Mi + s^Ma) + M|c2/3//i Q- Recall that E,i arises through the 
mismatch of soft terms between Lj and Hi as follows; 

The tree mass in Eq. p7|) gives the heavier mass scale, = Considering the 



atmospheric neutrino mass-squared difference, Am^^j^^ ^ 2.5 x 10~^ eV ~ m|, we get 



= 7.4 X 10-^ ^ — ^ (19) 
^ ^ VMz/ V0.05 eVy ^ ' 

Since the two mixing angles, 6*23 = 6'atm and 6*13, satisfy 

tan^23 = 6/^3 ~ 1 , I tan^isl = « 1 (20) 

we need ^1 < 0.3^2,3 to make small Q\z and ^2 ~ ^3 for near maximal atmospheric mixing. 
Thus, current neutrino oscillation data require 

\ (F^\^l^ ( ma ^1/2 



6 « 6 ~ 6 ~ 5.2 X 10-^- [-^\ — ^ (21) 
^ ^ CfiKMz) VO.OSeV/ ^ ' 



Let us now consider how one-loop corrections generate the neutrino masses and mixing 
accounting for the solar neutrino oscillation. In the bilinear model, the bi-large mixing of 
the atmospheric and solar neutrinos cannot be obtained under the assumption of universal 
soft terms 2- Thus, one needs to introduce non- universality in soft terms in order to 
accommodate the large solar mixing. 

Depending on the degrees of the deviation from the universality, we can consider two 
cases. First, the non-universality of soft parameters can arise due to small mismatches 
(likely to be caused by some threshold corrections) in the renormalization group evolution. 
In this case, the quantities Amf, Tn\.ij^/ei and fiABi are much smaller than the typical 
soft mass-squared so that the induced trilinear couplings in Eq. (jiup give the major 



contribution to the size of m2 ~ J Am: 



As discussed before, this causes the 



contradiction of Eqs. (0) and ((Tj). In other words, the condition of m2 ~ 8 meV yields 

ej ~ 4 X 10~^C/3 which is far above the first constraint in Eq. (fTBj) . 

However, we point out that there is a different way of reconciling the neutrino data with 

the baryogenesis requirement. Note that one cannot exclude the possibility that the loop 

mass is larger than the tree mass. For instance, one can take the superpotential bilinear 

parameter much larger than ^j, accepting a very small deviation of the non-universality 

or a cancellation among the terms see Eq. (fTHj) . In this situation, the heavier neutrino 

mass scale can be produced mainly by the bottom-sbottom loop which can be rewritten 

from Eqs. (jU and (fTU|) as follows: 

..loop _ 3/i| ml{Ab + fi tan (3) ^ 
Svr^ m( — mi mi 

As the above loop contribution determines the atmospheric neutrino mass and mixing, the 

condition (j2T| has to be replaced by 



..«.,«.3«8xl0-%^33^j . (23) 

Similarly to the previous discussions, €2,3 cannot satisfy the baryogenesis constraint ()16p at 
all. But, ei can be made arbitrarily small. Let us recall that it is sufficient to suppress 
lepton number violating couplings for just one lepton flavour. In our case, it is the electron 
number, which is implied by the smallness of 613. Now, in order for the tree mass (|17|) to 
produce the solar neutrino mass and mixing, we need 

1 /2 

6 ~ 6 ~ 3 X 10-^1 (^) . (24) 

C/3 \Mz/ 



For such small ei and the small deviation of the universality, we expect ~ 771 and thus 
the variables £11,12 in Eq. (fT^ can be approximated by 

neglecting M| terms. From this, one can see that the baryogenesis constraint can be 
satisfied ifl<tp< (m?^ — Tn\)/{ml — m^c^). 

Secondly, we consider the more general non-universality implying that Amf, m^L-H^/^i 
and fiABi are of the order m^. In this case, the neutral scalar and neutralino exchange 
loops can give important contributions to the one-loop mass as long as tan (3 is not too large 
and the large misalignment between C,i and rji is allowed. Adopting the result of Ref. 
the one-loop mass coming from the neutral scalar loops is roughly given by 



large 



Mj"^ ^ ^^x^^^^l'^M^lo^ml) (26) 

where Bo{x, y) = f ~ ^ + ^ ^^"^ represents the neural Higgs bosons, /i, H and 

A. Neglecting unimportant contribution of ^j, the variables 9is are approximately given by 



rjiS 13171^ 



2 rn^Ca^.p - M|c2^c»+^ 



{ml. - ml){m\ - mjj) 

{ml. -ml){m\ - mj^) 
2 



~ ir]iS(} 2" ^ 2 (27) 



mi - m^ 

,2 



where mh,H are the Higgs boson masses at tree- level determined by jj = 1/2 [m^ + M| =p 



'm\ + M|)2 — Am\M'^clp], and the angle a is defined by C2a = C2i3{m\ — M^) / {mj^ — m^j) 
and S2a = S2f3{m\ + Mz)/{m\ — m^j). Our convention for the pseudo-scalar Higgs boson 
mass is that m\ = —^B/cpSp. Requiring m2 ~ one obtains 

'300 GeVV^^ f f ^2 \ 



\ m^o J V"^x° / ^° meV/ 
As discussed in Ref. jl^, the large mixing of solar neutrinos require ~ 6*2,/,. This has to 
be contrasted the condition ^2 ~ ^3 ()21|) for the large atmospheric neutrino mixing. 

From Eqs. (fT^ . (PT|) and (|^. one sees that the couplings eji,j2 are required to be smaller 
than C,i or ^^^^ by one-order of magnitude. Thus, it is generally difficult to satisfy both the 



8 



baryogenesis constraints and obtain the realistic neutrino masses and mixing. However, 
it is not impossible to find some reasonable parameter space where both requirements are 
reconciled, which is due to the fact that the variables £:ii,i2 and C,i or have different 
dependencies on the input parameters. Comparing Eq. (jl4j) with Eq. (j27p . we notice that 
^ii,i2 (or and rji) can be made small while keeping ^ 6 x 10^^ (|28p when the sneutrino 
mass mj>^ is close to one of the Higgs boson masses, m/j, rriH and m^. Since the heavy Higgs 
scalar mass, mn, is usually close to the pseudo scalar mass, m^, and ~ 1, we find it 
better that the sneutrino mass is closer to the light Higgs scalar mass, that is, m,>j ~ m/^. 
Barring cancellation, both terms in en (fT^ should be less than 3 x 10^''. Again here, this 
is possible for i = I, that is, the electron number violating parameters, £11,12, can only be 
suppressed for our purpose. For illustration, let us calculate en + ^i, ei2 and for the 
cases with niA = 100, 300 GeV and tan/5 = 3, 30. In what follows, we present the values of 
i^ih, GiH, 6iA] + £i2) normalized with ei2 = 1, indicating the rough ranges of m,>. allowing 
for en > ei^/20: 

Case 1 tp = 3, mA = 100 GeV, ruh = 60 GeV (29) 

-77, +64, -9.1; -3.6, 1) for m^^ = 55 GeV 
-22, -21, -6.5; -0.96, 1) for m^^ = 71 GeV 
Case 2 = 30, tua = 100 GeV, ruh = 90 GeV (30) 

-113, +342, -77; -17, 1) for m^, = 73 GeV 
-66, -210, -49; +11, 1) for m^)^ = 115 GeV 
Case 3 = 3, = 300 GeV, rrih = 72 GeV (31) 

;+29, +53, -8.4; -1.5, 1) for m^^, = 60 GeV 
-10, -26, -5.5; +0.54, 1) for m^, = 90 GeV 
Case 4 = 30, tua = 300 GeV, ruh = 91 GeV (32) 

;+18, +289, -74; -14, 1) for m^^ = 75 GeV 
-6.9, -212, -49; +11, 1) for m^, = 115 GeV 

From the above calculation, one can see that the non-erasure condition can be satisfied if 
the difference between the sneutrino and the light Higgs boson mass is within 10%. In order 
to confirm the above properties, we made a numerical calculation to find a set of points 
satisfying both the baryogenesis constraints and the atmospheric and solar neutrino data. 



For this, we incorporate the exact formulae for the neutrino mass matrix derived in Ref . |15| . 
In Figures 1 and 2, we plot the variable eu in terms of the electron sneutrino mass m;>^ for 
all the points accommodating all the observed neutrino data for Cases 1 and 2. The plots 
clearly show the suppression of su when the sneutrino mass is close to a Higgs boson mass. 
Similar behavior is also found in Cases 3 and 4. 

Another way of suppressing en is to arrange a cancellation between two terms in en. 
From Eqs. ()14|) and ()27|) . one generally has 

£ii ~ -6 - t/3ei2 and ~ t/3ei2 (33) 

for mp- ^ Mz- Now, one can see that the conditions (fTB|) and (PHj) can be satisfied for 
t/3 ~ 30 with the cancellation in en ~ + 6i^. Again, this can work only for the electron 
direction with ~ since Eq. fETj) shows ^ ^^j,/, ~ 6 x 10~^ for i = 2,3 and large tan /3. 

In conclusion, we have investigated how the cosmological requirement for a successful 
baryogenesis can be reconciled with a realistic neutrino mass matrix in the R-parity vio- 
lating version of supersymmetric standard model. Our main focus has been to see whether 
the B — L violating interactions can be sufficiently suppressed in order not to erase a pre- 
existing baryon or lepton asymmetry of the universe. Such a baryogenesis constraint cannot 
be satisfied if the trilinear R-parity violating couplings are introduced to explain the atmo- 
spheric and solar neutrino masses and mixing under the assumption of the universal soft 
supersymmetry breaking masses. In the bilinear model, the observed neutrino data can be 
well explained if the non-universality is allowed. Our analysis shows that the non-erasure 
condition can be met by suppressing the electron number violating parameters, which is re- 
lated to the smallness of the angle 6'i3. In the case of a large violation of the universality, the 
electron sneutrino mass has to be nearly degenerate with the light Higgs scalar mass. For a 
small violation of the universality, we argued that the situation of the loop mass dominating 
over the tree mass is preferred contrary to the usual consideration. A consequence of our 
analysis is that the bilinear R-parity violating supersymmetric standard model can provide 
a framework not only for a realistic neutrino mass matrix but also for a successful baryoge- 
nesis through the Affieck-Dine Mechanism. Finally, let us note that our consideration is not 
relevant if the electro-weak baryogenesis is operative. 
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FIG. 1: The quantity en is shown as a function of the electron sneutrino mass for all the 
points generating the required neutrino masses and mixing for Case 1 in Eq. H29|) . 
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FIG. 2: Same as FIG. 1 for Case 2 in Eq. 
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